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1.  INTRODUCTION 

Adaptive  approaches  and  a-posteriori  estimates  of  the  accuracy  of 
computed  solutions  have  recently  attracted  great  Interest.  Nevertheless  the 
basic  notions ,  alms  and  goals,  principles  for  comparisons,  etc.  are  not  yet 
well  established.  We  try  In  this  paper  to  formulate  some  principle  notions 
and  approaches  related  to  this  new  field  of  numerical  mathematics  and 
computational  engineering  and  show  their  practical  importance.  We  will 
Introduce  the  notions  of  feedback,  adaptive  approaches,  a-posteriori 


estimates,  etc.  and  formulate  illustrative  theorems  In  a  model  setting. 
Although  the  available  theoretical  results  cover  a  much  larger  area  than 
formulated  In  this  paper,  the  general  theory  still  waits  to  be  developed. 

We  will  place  some  emphasis  on  the  analysis  of  the  one  dimensional  model 
problem  treated  by  the  h,  p  and  h-p  versions  of  the  flnltle  element 
method.  In  this  case  we  are  relatively  easily  able  to  formulate  the 
Illustrative  theorems  and  to  show  the  general  Ideas. 

The  results  mentioned  here  Illustrate  more  general  theory  and  Its 
numerical  Implementation  developed  by  a  group  at  the  University  of  Maryland. 

We  mention  In  this  connection  especially  contributions  of  W.  Gul  [1],  B.  Guo 
[2],  Ch.  Mesztenyl,  A.  Miller  and  M.  Vogellus. 

We  refer  also  to  [3],  [4],  [5],  [6],  [7]  for  some  relevant  approaches 
and  Ideas. 

Our  results  relate  mostly  to  linear  elliptic  problems.  For  nonlinear 
problems,  problems  of  other  type  such  as  parabolic  and  hyperbolic,  etc., 
some  other  aspects  need  to  be  emphasized;  nevertheless  we  will  not 
address  them.  For  some  of  these  aspects  related  to  the  solution  of  nonlinear 


problems,  see  [5],  [8]. 


THE  BASIC  MODEL  PROBLEM 


2.1,  The  one  dimensional  problem 

Let  us  be  Interested  In  the  one  dimensional  problem 


(2.1a) 

-(au')'  +  bu  *  f. 

x  €  I  -  (0,1) 

(2.1b) 

u(0)  -  u(l)  -  0 

where 

(2.1c) 

0  <  <Xq  <  a(x)  <  aj  <  «,  0  < 

b(x)  <  <Xj  <  •. 

We  will  assume  that  f  is  such  that  the  solution 

uq  of  the  problem  (2.1a,b) 

belongs  to 

H*(I)j  i.e.  uq  €  H*(I)  where 

(2.2)  H^I) 

•  {u(x)|lul?  *  /  (u'2-h»2)dx  <  •, 

0 

u(0)  -  u(l) 

"  0} . 

The  triple 

(a.b.f)  will  be  referred  to  as  input  data. 

We  denote 


(2.3) 


lul 


E(I) 


/  (au'2+bu2)dx. 
0 


The  solution  of  (2.1)  is  formulated  in  the  usual  weak  form:  Find  uq 

Oi 

H  (I)  so  that 


1  1 

(2.4)  B(un,v)  -  /  (au'v'+buv)dx  -  /  f  dx 

0  0 


holds  for  all  ▼  €  H*(I). 

A  special  consideration  will  be  given  to  the  case  when  the  solution  uq 


w  where 


(2.5) 


w(x) 


1*-y1 


[x-rl  +  CjX  +  C2, 
*  "  Y»  *  >  Y» 
0,  x  <  y. 


a  >  V2  not  integral  and  C^,  i  ■  1,2  are  chosen  so  that  (2.1.b)  holds. 

This  solution  (with  y  ■  0)  is  a  one  dimensional  model  of  the  solution  with  a 
singularity  caused  by  the  presence  of  cracks  and  corners  in  the  two  dimen** 
sional  setting. 

We  are  Interested  in  the  accuracy  of  the  approximate  solution  measured 
in  the  energy  norm  (2.3): 


<  tIu0IE  * 


2.2.  The  two  dimensional  problems 

Let  us  be  Interested  in  the  elasticity  problem  (plane  strain)  of  the 
cracked  panel  shown  in  Fig.  2.1.  We  assume  that  the  Poisson  ratio  v  ■  .3, 
Young's  modul  E  ■  1  and  the  tractions  are  such  that  the  solution  coincides 
with  the  solution  for  the  slit  half  plane  with  stress  intensity  factors  kj  - 

kII  “ 

The  energy  norm  (the  analog  to  (2.3))  and  the  weak  formulation 
(analogous  to  (2.4))  is  the  usual  one. 


Figure  2.1.  The  scheme  of  cracked  panel 
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The  second  problem  will  be  to  solve  on  S  ■  (x^^l  |xj|  <  j, 
I  *2 1  <  V2  }  the  boundary  value  problems 


(2.6a) 


Au^  ■  0  on 


(2.6b) 


on  8  S,  i  -  1,2 


Re[(z  -  -|)lg(z  -  -|)] 


Im[(z  -  4") lg( z  -  -5-)] , 


3.  THE  FINITE  ELEMENT  METHOD  (FEM)  FOR  THE  ONE  DIMENSIONAL  PROBLEM 


(3.1a) 


i  0  <  x$  <  xf  <  —  <  x£U)  -  1, 


(3.1b) 


m  \  wd  A  _  A 

vxj,xj+r*  nj  xj+i  xj» 


b(A)  - 


bj,  h(A)  ■  min  hj* 


(3.1c) 


PA  -  (pf »•••  »Pk(A)^»  PJ  >  l»  lnt«8Ml- 


A  la  called  meab,  x^  the  nodal  points,  Ij  the  elements  and  p^  the 
degrees.  We  denote  by  S(A,p  )  c  H1(I)  the  set  of  all  functions  which 
restricted  to  Ij  are  polynomials  of  degree  p^.  For  k(A)  •  (A,pA)  we 
define 


(3. Id) 


N(k(A))  -  N(A,p°) 


-  I  (p*+D  “  (M(A)+1). 

j-1  2 

N(k(A))  is  the  dimension  of  S(A,p^)  (i.e.  the  number  of  degrees  of 

freedom) . 

The  finite  element  solution 

u  ■  u(k(A))  "  u(A,pA)  €  S(A,pA) 

is  defined  In  the  usual  way: 

1  » 

(3.2)  B(u,v)  -  J  fv  dx  for  every  v  $  S(A,p  ). 

0 

The  formulation  of  the  finite  element  solution  Is  quite  general  and  Includes 
all  three  versions  of  the  method  namely  the  h,  and  h-p  versions  of  the 
finite  element  method. 

In  two  dimensions  the  formulation  Is  analogous  but  more  complex. 


4.  THE  FEEDBACK  FINITE  ELEMENT  METHOD.  THE  ACCEPTANCE  CRITERION 


By  the  feedback  finite  element  method  we  will  understand  a  specific 

^1 

construction  of  a  sequence  £  “  (k^,...),  ■  (A^,p  ),  and  the  sequence  of 


associated  finite  element  solutions  ■  u(k^)  1  -  1,2,...  . 

The  sequence  £  Is  usually  constructed  recursively  and  Is  called  a 
trajectory.  In  general  depends  on  all  u^,  A  -  1,2,...  1-1  and 

on  the  Input  data  (a,b,f). 

Us  will  write 


(4.1) 


cA  m  A(Kj,...,Kj_ji  Up...,u^_^,  a,b,f) 


and  call  A  a  transition  operator.  In  addition  we  will  assume  that 


jc_  e  Q,  i.e.  is  restricted  to  a  special  distributions  of  nodal  points  and 

degrees.  The  set  Q  is  given  a-priori  and  the  transition  operator  which 
leads  to  a  trajectory  £  €  Q  will  be  called  a  Q-transition  operator. 

To  every  an  error  estimator 

(4.2)  ■  E  ^(tc^.u^  »a»b  »f ) 
is  defined  and  the  solution  u^  is  accepted  if 

(4.3)  E  i(<i,u1,a,b,f )  <  t»«ilE(I) 

where  t  is  an  a-priori  given  tolerance. 


5.  PRINCIPLES  OF  THE  SELECTION  OF  A  TRANSITION  OPERATOR 

Obviously  there  are  many  possible  transition  operators  which  could  be 
considered  in  practice.  An  optimal  selection  depends  on  many  factors,  th-i 
class  of  problems,  computer  technology,  etc.  Nevertheless,  it  is  of  ma jc  - 
interest  o  introduce  notions  which  will  allow  one  to  characterize  basic 
properties  of  feedback  procedures,  to  make  comparisons  and  to  develop  a 
systematic  mathematical  theory. 

5.1.  Adaptivity  with  respect  to  the  convergence  measure 

Given  the  transition  operator  A,  the  problem  is  to  describe  the  set 
Rj  of  the  input  data  (a,b,f)  and  the  Initial  meshes  so  that 

(5.1)  ,uo“'j(<i),g(i)  0  as  i  ♦  •. 

Such  an  operator  A  (and  the  feedback  finite  element  method)  will  be  called 
adaptive  relative  to  ^  with  respect  to  the  convergence  measure.  Obviously, 
in  comparing  two  transition  operators  from  this  point  of  view,  we  prefer  the 


one  Whose  set  Rj  is  larger  or  covers  the  set  of  input  data  which  are  more 
relevant  In  applications.  Any  feedback  finite  element  method  used  in  practice 
has  to  be  adaptive  with  respect  to  the  convergence  measure  for  a  large  set 
Rj.  Therefore,  a  mathematical  characterization  of  Rj  is  of  major 
importance. 

5.2.  Adaptivity  with  respect  to  the  convergence  rate  measure 
Define 


(5.2) 


®(N)  -  ®(N,Q,a,b,f) 


inf 

tc€Q 

N(k)»H 


»u0-S(<)l 


E(I) 


and  call  $(N)  the  rate  function.  Obviously  ®(N)  is  the  smallest  possible 

error  obtainable  by  the  mesh  of  any  trajectory  of  •»  with  the  degree  of 
freedom  N.  (We  remark  that  $(N)  has  not  necessa  'ily  been  defined  for  all 
N.) 

Given  the  Q-transltlon  operator  &  and  a  number  p  >  1,  the  problem 
is  to  describe  the  set  R2(p)  of  the  input  data  (a,b,f)  and  the  initial 
meshes  ic^  so  that 


(5.3) 


sup 

1-1,2,... 


l"(lCi)-U0,E(I)  . 
®(n(k1)J -  < 


P 


or  the  set 


R 


2 


(5.4) 


U  R,(p) 


The  analysis  of  the  rate  function  <$(N)  for  various  sets  of  input  data  and 
the  admissible  set  of  trajectories  jc  €  Q  is  of  major  importance.  It  gives 
a  guide  for  the  design  of  an  effective  transition  operator. 


The  Q-transition  operator  A.  (reap,  the  feedback  finite  element 
method)  satisfying  (5.3)  (resp.  (5.4))  is  called  p-adaptive  (resp.  adaptive) 
relative  to  R?(p)  (resp.  fy)  with  respect  to  the  convergence  rate 
measure. 

Obviously  the  definition  is  strongly  tied  to  (5.3),  where  the  error  of 
the  ideal  mesh  is  in  the  denominator.  We  can  weaken  (5.3)  in  order  to 
increase  and  still  practically  require  that  the  performance  of  { 

to  be  nearly  ideal.  For  example,  we  can  ask  instead  (5.3)  that 


(5.5) 


,3<*i>“VK 

8  ip  - - 

i-1,2,...  [«(N(ie  ))]P 


<  P 


with  0  *  1  -  e  and  e  >  0  : ..pall* 


5.3.  Adaptivity  with  respect  o  the  work  measure 

Let  us  define  the  computational  effort  of  determining  the  finite  element 

solution  u(k1)  by  an  expression  W(k^).  We  can,  for  example,  take 
0 

W(<^)  ■  [N(k^)]  where  0  is  a  properly  chosen  number. 

MS  MS 

Computation  of  u(tc^)  assumes  that  u(kj),  j  *  l,...,i-l  were  already 
determined.  Therefore,  the  total  work  W(k^)  needed  in  determining  ufc^) 
is: 

_  i 

(5.6)  W(ic<)  -  l  W(k,). 

j-1  2 

We  Introduce  now  the  proportionality  factor  x(*^): 


I 


The  problem  now  is  to  describe  the  set  ^(p)  of  input  data  and  the  initial 


meshes  so  that 


(5.8) 


sup  x(0  <  P* 

i-1,2 . 


The  transition  operator  A  resp.  the  feedback  finite  element  method 
satisfying  (5.8)  is  called  p-  adaptive  with  respect  to  the  work  measure 
relative  to  £3(0). 

5.4.  The  adaotivitv  with  respect  to  the  effectivitv  index  of  the  error 


estJ  .ator 


(5.9) 


9(Xj) 


Ei(lCi,“i,a,b,f) 

,u0_“(lfi),E 


©(x^)  is  called  the  effectivity  index  of  the  estimator.  Given  a  transition 
operator  A,  the  problem  now  is  to  describe  the  set  R^p)  of  the  input  data 
and  the  initial  meshes  Xj  such  that 


(5.10) 


f  <  00^)  <  p 


and  the  set  R4  such  that 


(5.11] 


!©(x^)-l|  +0  as  i  -*•  ®. 


If  in  addition  ©(x^)  >  1  then  we  will  say  that  R^  is  an  upper 
estimator.  The  problem  now  is  to  describe  the  set  ^(p)  (resp.  R4).  The 


The  notions  of  adaptivity  for  two  dimensional  problems  are  defined  in 


analogy  to  those  for  one  dimension. 

The  notions  of  feedback  and  adaptivity  similar  to  as  introduced  above 
were  first  suggested  in  [7]. 

5.6.  The  selection  of  the  feedback  finite  element  method 

We  have  introduced  some  basic  notions  which  allow  one  to  compare  various 
feedback  FEMs  and  error  estimates.  Of  course  other  important  comparison  could 
be  considered.  Nevertheless,  we  believe  the  goals  of  the  feedback  procedures 
always  have  to  be  exactly  defined  and  theoretically  and  experimentally 
investigated. 

6.  RATE  FUNCTION  9  AND  ITS  PRACTICAL  MEANING 

The  rate  function  9  introduced  in  Section  5.2  is  of  principal 
importance  because  it  describes  the  performance  of  an  ,,ldeal,,  feedback 
procedure.  It  gives  also  guidelines  for  the  design  of  concrete  feedback  FEM 
by  allowing  one  to  judge  consequences  of  programming  and  other  considerations 
(influencing  e.g.  the  set  Q). 

6.1.  One  dimensional  model  problem 

We  now  mention  here  some  typical  illustrative  results  in  the  case  when 

ug  ■  w  with  w  given  by  (2.5). 

Theorem  6.1.  Let  y  <  1  in  (2.5)  and  a,  b  satisfy  (2.1c).  Let  Qj  be  the 

A  v  A 

set  of  all  meshes  such  that  h(A)/h(A)  <  a  <  •  and  pa  *  (p,*».,p),  p  >  1 


•*.  «*. 


(l.e,  consider  the  classical  h-version  with  quasiuniform  mesh  and  elements  of 
degree  p).  Then 


(6.1)  <  C2(a,Y,p,ff)N~p 

where 

a)  for  0  <  y  <  1: 

(6.2)  P  “  min(p,a-j), 

b)  for  y  ^ 

(6.3)  p  ■  p. 


Theorem  6.2.  Let  0  <  y  <  1  and  a,  b,  satisfy  (2.1c).  Let  Q2  be  the 
set  of  all  meshes  and  pA  ■  (p,...,p),  p  >  1  (l.e.  we  consider  the 
classical  h-version  with  arbitrary  meshes  and  elements  of  degree  p).  Then 

(6.4)  C1(a.Y.p)N_P  <  tfN,^)  <  C2(a,Y  ,p)N-p. 


We  see  an  essential  downgrading  of  the  performance  of  the  h  version  if 
we  restrict  ourselves  to  the  quaslunlform  meshes  and  the  solution  is  not 
smooth.  In  the  case  of  arbitrary  meshes  the  performance  is  governed  here  only 
by  the  degree  p  and  not  the  smoothness  of  the  solution  (characterized  by  a). 

(6.4)  shows  also  that  use  of  higher  degrees  is  advantageous.  (We  will 
see  this  further  in  connection  with  the  p- vers ion  and  by  the  numerical 
examples  in  Section  7.3.)  Nevertheless,  we  have  to  assume  here  that  the 
achieved  accuracy  is  sufficiently  good,  l.e.  we  are  in  an  asymptotic  range. 
This  could  be  practically  impossible  if  a  is  close  to  ®  ■  *55 

and  we  are  restricted  by  double  precision  to  implement  proper  meshes  without 
round-off  problems  as  will  be  seen  in  Section  7.3). 
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Theorem  6.3.  Let  y  *  0  and  a,  b  satisfy  (2.1c).  Let  Q3  be  the  set  of 
all  meshes  and  all  distribution  of  p  (i.e.  we  consider  the  general  h-p 
version).  Then 


(6.5) 


C^a)  /(a  -  V2  )N 


<  ®(N,Q3)  <  C2(a)r 


/  («  -  V2 )N 


where 


(6.6) 


-  (/T  -  l)2,  0  <  t  <  |  a. 


We  see  that  the  h-p  version  leads  to  an  exponential  rate  with  respect 
to  /N. 

Let  us  address  now  the  p-version.  To  show  essential  featuraes  we  will 
restrict  ourself  to  the  case  when  only  one  element  is  used,  i.e.  M(A)  ■  1. 

Theorem  6.4.  Let  a  ■  1,  b  ■  0  and  <J^  be  the  set  of  all  meshes  with  one 
element.  Then 

a)  for  y  *  0: 

(6.7)  ®(N,Q4) 

b)  for  y  <  0: 

Let 

q  _  /1-y  -  /— y 

/1-y  +  /-y 


T(g;  Isin  x(g-l) |  1 


N 


2g-l 


(i+o(l)) 


then 


bl)  for  0  <  q2  <  1  - 


14 


(6.8) 


C.  N-a  .  ,  ct-  V2 

q  (1-q  >  }  <  WO*) 


/  1-q2  N  N®”^2 


<  C,(l  +  — i— )  £1  {_!_  +  ( l-,2)*'  l/2  ) 


✓  1-q' 


N®  W2 


VT 


b2)  for  1  -  4  <  q2  <  Is 


N 


(6.9)  C  t  {~+  (l-q2)®"1/2  }  <  *(N,Q, ) 

V1/2n«-1/2 


N-a 


<  C2^— d-q2)a‘1/2} 
Na-V2Na-l/2 


where  Cj  and  C2  are  positive  and  independent  of  N  and  q 
c)  for  0  <  y  <  Is 


(6.10) 


C(a,Y>  — -1-—  <  ®(N,Q4> 

Na-V2 


N 


2a- 1 


for  0  <  6  <  — 

N 


<  C(a,y) 


f  l)  a-  V2 

'•N; 


for  6  <f 


where  6  ■  min(0,x-0),  0  ■  arccos  (2y-l).  For  y  <  0  the  rate  is 

exponential,  for  y  ■  0  the  rate  is  twice  than  that  for  0  <  y  <  1. 

Comparing  ^(NjQ^)  and  $(N,Q^)  we  see  that  the  performance  of  the  p- 
version  is  never  worse  than  the  h-version  with  arbitrary  p  for  the  uniform 
mesh  but  it  can  be  twice  (y  *  0)  or  more  (y  <  0)  as  good.  This  is  true 


SS 


►.V-. 


r.  *  ,  • 

>v  A 

■N*  V* 


(6.13) 


r  .  /  Ig1^/  lg</2  -1) 
r*  r 

where  r  Is  given  by  (6.6). 

Theorems  of  the  mentioned  type  are  Important  for  the  design  of  feedback 
procedures.  For  example,  although  the  estimate  (6.11)  is  conservative,  we  see 
that  if  the  optimal  positions  of  the  nodal  points  are  not  obtained  by  at  most 
one  or  two  iterations,  it  is  more  effective  to  increase  the  number  of  nodal 
points  with  the  use  of  recursive  meshes  than  to  further  seek  the  optimal 
position. 

6.2.  The  two  dimensional  model  problem 

The  situation  in  two  dimensions  is  similar.  As ' Illustration  we  will 
mention  theorems  which  deal  with  our  first  model  of  elasticity  problem. 

Theorem  6.7.  Consider  the  h-version  with  quasi-uniform  mesh.  Then 

(6.14)  C^p^N-1^  <  *(N)  <  C2(p,o)N~  ^  # 

Theorem  6.8.  Consider  the  h-version  with  arbitrary  meshes.  Then 

(6.15)  C1(p)N“P/2  <  ®(N)  <  C2(p)H"p/2  . 

Ue  see  that  Theorems  6.7  and  6.8  are  analogous  to  6.1  and  6.2  where  N 
is  replaced  by  fa  . 

As  the  analog  of  Theorem  6.3  we  have 

Theorem  6.9.  Consider  the  hrp  version,  Then 

3_ 

®(N)  <  C  e”p/N  , 


(6.16) 


6  >  0 


There  Is  no  sharp  estimate  for  g  available,  neither  Is  It  known  whether  the 
cube  root  In  (6.16)  Is  the  best  possible. 


In  contrast  to  the  one  dimensional  case,  a  restriction  can  be  placed 
here  not  only  on  the  meshes  but  also  on  the  structure  of  polynomials  of  the 
elements.  It  can  be  shown  that  g  In  (6.16)  can  be  Increased  If  the 
polynomial  elements  do  not  necessarily  consist  of  all  shape  functions  of  the 
same  degree  (l.e.  are  chosen  selectively)  in  comparison  with  the  case  when  all 
shape  functions  of  the  same  degree  are  used.  This  fact  could  be  Important  for 
the  design  of  a  feedback  procedure.  See  e.g.  [9]. 


7.  SOME  SPECIFIC  FEEDBACK  PROCEDURES 

We  will  formulate  a  specific  feedback  finite  element  method  and  analyze 
its  properties  In  the  light  of  Section  6. 


7.1.  The  feedback  procedure  for  the  h-verslon 

Let  us  consider  the  problem  (2.1a,b,c)  and  for  simplicity  use  only  p  ■ 
1.  Let 


(7.1) 


(hj) 


(az'^+  bz^)dx 


J-l 


where  z(x)  is  such  that 


(7.2)  -*(az')'  +  bz  -  f  -  [-(a Z'Y  +  bS]  -  r(x) 

(7.2a)  z(Xj-i)  "  "  °’ 

z(x)  is  the  solution  of  the  original  problem  on  Ij  when  f(x)  is  replaced 
by  r(x).  We  will  assume  that  the  number  rf j  can  be  well  approximated  by  rjj 
which  is  easily  computable.  For  example  we  define 


(7.3) 


(*? 


i - 7  J  r  C*)dx 


12  a(xA  )  A 

j-  W  *3-1 


j”  lk 


xA  +xA 

3-13 


The  number  Tfj  Is  called  an  error  indicator.  The  error  estimator  E  la  now 
defined  as  follows: 


(7.4) 


r  r^,2- 

3-1  j 


We  will  assume  that 

lhj-^1  <  4»(h j)  Ihjl 

where  4>(hA)  >  0  as  hA  -*■  0  which  allows  us  consider  T) j  and  as 

identical. 

A, 

***  i 

Assume  now  that  A^,  u^,  T)  ^  have  been  computed.  Then  A^+1  Is 

Ai  Ai 

construced  by  bisecting  all  elements  of  A^  for  which  hj  >  X  max 

where  0  <  \  <  1  is  an  a-priori  given  parameter  which  steers  the  process. 
We  will  now  state  basic  theorems  related  to  this  feedback  procedure. 

Theorem  7.1.  The  set  R|  (see  Section  5)  consists  of  all  input  data 

o, 

satisfying  (2.1c),  all  solutions  of  H  (I)  and  an  arbitrary  initial  mesh 
Theorem  7,2.  For  any  0  <  X  <  1,  ^  c  ^i>  ^2  **  ^1* 


Theorem  7.2  shows  that  to  get  adaptivity  with  respect  to  the  convergence 
rate,  stringent  restriction  on  R2  has  to  be  made.  In  [4]  we  formulated 
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sufficient  conditions  so  that  (a,b,f)  (  ty. 

As  an  illustration  we  state 

Theorem  7.3.  Let  (2.abc)  be  satisfied*  Assume  in  addition  that  a(x) 
satisfies  Holder  condition  with  exponent  0  >  y.  Then  if  0  <  X  <  1,  the 
input  data  (a,b,f)  such  that  uq  -  v  and  arbitrary  initial  mesh  Aj  belong 
to  i?2* 


Combining  Theorems  7.2,  7.3  and  6.2  we  see  that 

luQ-i^lg  <  CCy.oOn”1  . 

Under  what  condition  the  input  data  belong  to  ^(p)  with  an  a-priori 
given  p  is  not  known. 

Theorem  6.1  shows  immediately  that  if  X  -  0,  y  >  0  and  a  <  3/2 

mm 

then  uq  ■  w  does  not  belong  to  Nevertheless,  if  y  <  0  end 

0  <  X  <  1,  uq  ■  w  belongs  to 

Theorem  7.4.  Let  X  ■  1.  If  uq  ■  w  and  0  <  y  <  1,  then  the  input  data  do 
not  belong  to  Rj(p)  for  any  p.  If  0  <  X  <  Xq  <  1  and  a,  b  satisfy  the 
assumptions  of  Theorem  7.3,  then  uq  ■  w  belongs  to  Rj(p)  with  p  <  pQ(a,y). 


Theorem  7.5.  Let  conditions  (2.abc)  hold,  a' ,  b  be  continuous  and  f 
L2(I).  Then 

10(0-11  <  Ch2(A) 

with  C  dependent  on  the  input  data  but  Independent  of  the  mesh. 

Any  design  of  a  feedback  approach  has  to  involve  a  compromise  among 
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various  requirements.  Theorems  of  the  type  mentioned  above  (in  larger 
generality)  proven  theoretically  or  experimentally  (as  conjectures)  allow  one 
to  make  a  rational  choice  and  an  assessment  of  the  design. 


7.2.  The  feedback  procedure  for  the  h-p  version 

Analogously,  but  in  a  much  more  complex  way  a  feedback  procedure  for 
the  h-p  version  can  be  designed.  The  main  difficulty  is  to  decide  when  to 
refine  the  mesh  and  when  to  increase  the  degree.  We  cannot  go  into  details 
here,  but  remark  that  we  have  designed  a  feedback  approach  and  analyzed  its 
properties.  The  results  will  be  shown  in  numerical  examples. 


7.3.  Numerical  examples  for  the  one  dimensional  problem 

In  this  section  we  will  show  numerical  data  and  analyze  them  from  the 
theoretical  point  of  view.  He  will  also  mention  in  this  connection  additional 
theoretical  results  pertinent  to  the  analysis. 

Let  e  ■  uq  -  Q  and 


lei 


R 


lei, 


,U0*E 


and  consider  the  performance  of  the  finite  element  method  when  Uq  ■  w  and 
w  be  defined  by  (2.5).  We  consider  a  ■  1,  b  ■  0  in  (2.1). 

Riga*  7.1  ab  show  the  performance  of  the  h-verslon  with  uniform  mesh 
for  y  •  .3  and  c  •  .7  and  1.1  (in  log  log  scale).  The  shown  slope  is 
the  theoretical  one  based  on  Theorem  6.1  (p  ■  .2  and  p  -  .6). 

Figaros  7.2  ab  depict  the  performance  of  the  feedback  h-verslon  (for  p 
•  1,2)  and  the  feedback  h-p  version. 


We  see  that  the  rate  of  convergence  of  the  feedback  h-verslon  is  0(N-i) , 
and  0(N  )  for  p  »  1  and  p  ■  2  independently  of  a.  Recalling  Theorems 

6.2  and  7.3  we  see  that  the  method  is  adaptive  with  respect  to  the  convergence 
rate  measure.  The  feedback  h-p  version  shows  exponential  rate. 

Table  7.1  shows  the  coordinates  of  the  nodal  points  constructed  by  the 
feedback  h-verslon  (p  *  1)  and  M  ■  80.  We  see  that  the  points  are  located 
very  close  to  the  singularity  point  (.3)  which  prevents  to  get  higher 
accuracy  due  to  limited  computer  accuracy. 


TABLE  7.1 


The  locations  of  the  nodal  points  for  the  feedback  h-version 
(p  *  1,  M  “  80,  y  *  .3,  a  *  .7) 


n 

X*  1 

*n 

— 

n 

yA 

*n 

0 

0.00  ! 

36 

.3000183105468750 

1 

.25  | 

• 

2 

.28125  ! 

46 

.301025390625 

3 

.286875  { 

47 

.3046875 

4 

.298828125  ; 

• 

5 

.2998046875  i 

66 

.390625 

6 

.2999267578125 

• 

7 

.2999877929875 

71 

.50 

8 

.29999542236283 

72 

.53125 

9 

.2999992370605468 

• 

• 

80 

.93750 

19 

.2999999999956344 

81 

1.00 

20 

• 

.3000000000029104 

23 


Fig.  7.4  shows  in  seal logarithmic  scale  (loglelR  and  /N)  the 

performance  of  the  h-p  feedback  FEM  for  a  ■  .7  and  a  -  1.1.  Using  Theorem 

6.3  (in  slightly  generalized  form)  the  graphs  in  the  asymptotic  range  should 
be  linear.  The  slopes  b,d  and  ac  shown  in  the  figure  are  based  on  (6.5)  and 
(6.12).  We  see  very  good  agreement  from  which  we  conclude  that  the  feedback 
is  adaptive  with  respect  to  the  convergence  rate. 

Fig.  7.5  is  illustration  of  Theorem  6.4.  It  depicts  the  performance  of 
the  p-version  (a)  with  one  element  only  for  various  y.  We  see  large 
sensitivity  to  y  of  the  accuracy  for  higher  p.  The  figure  also  shows  the 

performance  of  the  h  resp  h-p  feedback  FEM  (graph  b  resp.  graph  c) 


Y  “  .3,  «  -  .7,  p  -  1 


Y  ■  .3,  a  ■  .7,  p  ■  2 


Figure  7.3.  The  performance  of  the  h-verson  with  the  optimal  mesh  and  the 
feedback  h-version.  a:  the  feedback  h-verslon,  b:  the  h-version  with 
optimal  mesh. 


Figure  7.4.  The  performance  of  the  h— p  feedback,  b,d:  slopes  based  on 
(6.5),  a,c:  slopes  based  on  (6.12). 


Figure  7.5.  The  performance  of  the  p  and  feedback  h  and  h-p  versions 
for  y  <  0*  a!  p-version,  b:  feedback  h-version,  c:  feedback  h-p 
version. 


7.4.  Feedback  procedure  for  the  h-version  in  two  dimensions 

During  the  last  few  years  a  feedback  program  FEARS  (with  p  ■  1)  was 
developed  at  the  University  of  Maryland.  For  more  about  FEARS,  see  e.g.  [10] 


FEARS  has,  among  others,  the  following  properties: 


a)  It  has  been  proven  that  it  is  adaptive  with  respect  to  the  conver¬ 
gence  measure  for  a  wide  class  of  input  data* 

b)  It  is  conjectured  that  it  is  adaptive  with  respect  to  the  convergence 
rate  measure  for  a  large  class  of  input  data  of  practical  impor¬ 
tance.  Although  there  is  much  experimental  evidence  to  support  this 
conjecture,  the  proof  is  not  available. 

c)  Experience  shows  that  FEARS  is  adaptive  with  respect  to  the  work 
measure  with  p  <  2  or  3  in  all  practical  cases.  The  theoretical 
proof  is  not  available. 

d)  It  has  been  proven  that  FEARS  is  adaptive  with  respect  to  the 
effectlvlty  index  of  the  error  estimator  for  a  p  >  1  with  p 
independent  of  the  mesh  and  domain,  and  large  class  of  input  data. 

e)  It  has  been  proven  that  GOc^)  ♦  l  under  certain  relatively 
restrictive  assumptions.  There  is  a  conjecture  that  these 
assumptions  are  not  necessary. 

f)  It  is  conjectured  that  the  estimator  is  an  upper  one  for  a  large 
class  of  input  data  provided  that  the  achieved  accuracy  is 
sufficiently  good.  The  mathematical  proof  is  not  available,  but 
experimental  results  nevertheless  support  this  conjecture. 

7.5.  Numerical  Illustration 

We  solved  by  Fears  the  cracked  panel  elasticity  problem  introduced  in 
Section  2.2.  Table  7.2  shows  the  sequence  of  meshes,  constructed  by  FEARS, 
the  true  (relative)  error  and  the  effectlvlty  index.  Fig.  7.6  shows  the  mesh 


with  N  -  617 


TABLE  7.2 


The  performance  of  FEARS  computation. 


MESH 

DOF 

lei- 

“  r 

0 

"Ve  ' 

1 

101 

26.38 

.885 

2 

143 

21.35 

.991 

3 

221 

16.79 

1.058 

4 

301 

13.61 

1.116 

5 

617 

9.63 

1.088 

Figure  7.6.  The  mesh  constructed  by  FEARS  for  N  -  617. 

Table  7.2  shows: 

a)  The  rate  of  convergence  Is  0(N“  ^ )  and  supports  the  conjecturae 


that  FEARS  is  adaptive  with  respect  to  the  convergence  rate  measure 


b)  The  error  estimator  Is  an  upper  one  when  the  error  is  sufficiently 
small. 

c)  FEARS  is  adaptive  with  respect  to  the  effectivity  index  of  the 
estimator. 

d)  Because  /N^  ■  1.3  to  2,  and  it  has  been  seen  that  the  work 
(machine  time)  is  proportional  to  N1*8,  FEARS  is  adaptive  with 
respect  to  the  work  measure  with  p  “2  in  our  particular  case. 

8.  THE  MULTIPLE  LOAD  PROBLEM 

So  far  we  addressed  only  the  problem  of  adaptive  procedures  for  one  load 
only.  In  practice  we  have  often  consider  few  loads  (boundary  conditions)  and 
for  obvious  computational  reasons,  we  wish  to  use  one  mesh  for  all  these 
cases.  Theoretically  (and  practically)  the  multiple  load  problem  does  not 
bring  any  difficulties.  Vfe  will  show  it  in  a  second  example  in  Section  2.2. 

We  (theoretically)  solve  a  (noncoupled)  system 

du, 

■  °*  8T  *  ‘l 

du« 

■  °>  W  •  *2 

and  define  (for  example) 

-  i.,i|  +  i^i|  . 

(Of  course  we  determine  uj,  by  using  one  stiff nes  matrix  only).  Figure 
8.1a  shows  the  mesh  (constructed  by  FEARS)  for  (N  -  538)  (the  mesh  for 
U2  is  the  same  but  rotates  by  90°).  Figure  8.1b  shows  the  mesh  for  multiple 
loads  g},  g£,  N  "  531.  Figure  8.2  shows  le^lg  for  the  mesh  constructed 
for  Uj  as  simple  load  case  1  ■  1,2  end  for  the  mesh  which  is  used 
simultaneously  for  both  cases. 


Fig.  8.2.  The  error  I  e.  I E  using  single  load  and  multiple  load  meshes 
a:  single  load;  b:  multiple  load. 

9.  FEEDBACK  PROCEDURES  WITH  RESPECT  TO  OTHER  ACCURACY  REQUIREMENTS 

In  the  previous  section  we  concentrated  exclusively  on  the  accuracy 
measured  in  the  energy  norm.  Nevertheless,  in  structural  mechanics  we  are 
often  interested  in  the  computation  of  some  critical  values  (functionals)  as 
for  example,  stresses  in  some  points,  stress  Intensity  factors,  etc. 

There  are  various  ways  to  extract  these  value  from  the  finite  element 
solution.  These  extractions  are  implemented  by  the  postprocessing 
technique.  The  adaptivity  and  the  a-posteriori  estimates  for  the 
postprocessing  approach  have  some  special  features.  We  cannot  go  here  more  in 
details  and  we  refer  only  to  [3]  [II]. 
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